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unstable orbits. 
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It was observed in [1] that certain orbits for particles in a Schwarzschild space-time were 



O 

^\ 'unstable' to perturbations. A slight change in the initial conditions or perturbations of 

these orbits have exponential growth, in other words, one can find Lyapunov exponents. 
We study these orbits and the modifications of the Lyapunov exponents under effects of 
qq \ quantum gravity. 

O ' We find semiclassical corrections to the metric using coherent states in Loop Quantum 



Gravity (LQG), described in [2-5]. The corrected metric is then used to find unstable orbits. 
The Lyapunov exponents for circular orbits in the corrected metric are determined using 
the method of [1]. It is found that non-zero Lyapunov exponents occur even for normally 
classically stable orbits, but these cause insignificant effects for macroscopic astrophysical 
black holes. The time scale for the instabilities becomes significant (~ 1 second (s)) for 
primordial black holes with horizon radius of the order of 10 14 Planck length (l p ). It is 
then shown that this instability will not affect 'gravity waves' from these primordial black 
holes. However, this 'new' semiclassical instability might affect black hole mergers and 
other phenomena like quasinormal ring down from these primordial black holes [6, 7]. Next, 
we find the growth of the semiclassical corrections to classical circular orbit radii as per 
the classical Lyapunov exponents. This effect is finite for astrophysical black holes. The 
deviations grow to the order of 1 cm in a time of 1 second for black holes with horizon radius 
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10 42 / p which thus includes astrophysical black holes. However, in an actual accreting disc it 
will be difficult to isolate these instabilities from external matter perturbations. 

In the next section semiclassical corrections to the metric using coherent states of LQG 
are described. Section III finds the effects of these corrections on the unstable orbits. The 
last section describes some physical systems where one can detect this effect of quantum 
gravity. 



We discuss semiclassical states which are special as expectation values of operators are 
closest to their classical values in these states. In case of quantum gravity such special states 
were identified in LQG framework in [2, 3] using SU(2) coherent states defined by Hall. 

LQG is a formulation rooted in canonical gravity, where one takes a suitable time slicing 
of a space-time by identifying a fiducial time like direction. The induced metric in the 
time slices is denoted by q a j, (a,b=1..3) and the extrinsic curvature with which this slice 
is embedded in four dimensions is denoted by K a j,. The theory of loop quantum gravity 
is defined in terms of the triads or the square root of the metric q a b = e l a e\ and the 'spin 
connection' = e ,K ejV a ^bK ■ I — 1, 2, 3 represents an internal SU(2) index. A redefinition 
of the variables in terms of tangent space densitised triads E T a and a corresponding gauge 
connection A ! a simplifies the quantisation considerably. 



(e„ are the usual triads, the associated spin connection, and /3 the one parameter ambi- 
guity which remains named as the Immirzi parameter and is set to 1 for the discussions of 
the paper.) 

The 'holonomy' h e (A) = V exp(f e A.dx) and the momentum = f Se *E^ as integrals 
over an edge e and 2-surface S e (which is intersected by the edge e at one point) are defined 
as 'discretisations' of the gauge connection and the triads. The discrete variables give a well 
defined Poisson bracket. One can find a kinematic Hilbert space which carries a representa- 
tion of this algebra. A semiclassical state ij^ defined in the kinematic Hilbert space is given 
as a function of h e , g e = e lTlp ' h e (T 1 are the three 2x2 SU(2) generator matrices) and t (a 
semiclassical parameter). For the purposes of this discussion the semiclassical parameter is 



II. SEMICLASSICAL CORRECTIONS 
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taken as t = -ft = 10 , where r g = 10 n l p is the radius of the horizon of the Schwarzschild 
black hole, t — > represents the 'classical limit' and the wavefunction is given by 

1>* = E (2J + l)e" M XjW (2) 

3 

(j— 0,| ,1.., and Xj is the character of the jth SU(2) representation of the matrix gji' 1 ). 
The wavefunction is nicely peaked at the classical values of h e and P/. In the particular 
coordinates used to describe the embedding of the edges in the classical metric, they are 
labelled by e a , where a labels the three space dimensions. The fluctuations over the classical 
values can be obtained as a series in powers of t. ( We use the t notation to avoid confusing 
this with the time coordinate.) Thus 



<PJ>= < ^ > =PJ 



tf(P) (3) 



11^ 

where f(P) is a function and measures the first order correction to the classical value. The 
details of the function f(P) have been calculated previously in [3, 4]. Here we redo the 
calculations to ascertain the exact nature of the corrections. The action of the operator 
on the coherent state is given by 

p ^4(tXJ^ %) (4) 

(7 being a real parameter). The numerator of the expectation value of the momentum 
operator P/ 

where cosh(2;) = |Tr(e _7T/ gg T ), j = 0,1/2,1 labels the SU(2) angular momentum. By 
defining n — 2j + 1, one can rewrite the sum as 

In the limit t — > the above sum is converted into an Riemannian integral 



2i (~7~) 6 1, 1 I xe ^ sinh(za:)da; (7) 



^7/ 7=0 sinh(2;) Jo 
,c?7/ 7=0 4sinh(2) 



2 ^(4~) T^TT^ x pf?) ( 8 ) 



i 



where z = 2z/vt. Using the approximations of [3], one writes z = P Ca + 5 where P ea = 
PI PI and 5 = — ^TiiT 1 gg T ) / (2 sinh(P ea )) to first order in 7, and thus one can easily take 

: 2^^/ 4 e^A and 

(9) 



the derivative wrt 7. The norm of the wave function is H^ 4 



Pi 



I sinh P e 



Thus after some algebra 



< Pi >= Pi 



i + w br-coth(P ea ; 



(10) 



Interestingly in this approximation, the correction is proportional to the original value of 
Pi . The corrections are negative or positive as per the sign of Pi . Thus as per the above 



/(P) = p-(p--ooth(P ea ) 



(11) 



The sign of which is negative for all values of P £a 

In [4], I used a coherent state defined on a flat slicing of the Schwarzschild metric, and 
one starts with a metric defined in Lemaitre coordinates: 

dR 2 r3 i 4 /3 



ds 1 



-dr z + 



2/3 



+ 



r 2 J 3 (d9 2 + sin 2 



(12) 



The + sign in the above metric corresponds to a white hole with signals emitted from the 
centre r = 0, but — sign corresponds to a black hole with material falling to r = 0. Of 
course in a collapsing situation, the — sign should be solely relevant. The R, r coordinates 
are related to the Schwarzschild coordinates t, r using the following transformations: 

1 



l—dr = (dR±dT) dt = - 
r n 1 



-(dr±fdR) /'=-£ 



(13) 



The constant r = r c surface then has a flat metric. For the either of the signs used in the 
Lemaitre metric, the induced metric in the r = constant slice is ds 2 = dr 2 +r 2 (d8 2 +sin 2 8d<f) 2 ) 
where r l l 2 r~ l l 2 dr = dR is the transformed coordinate on the slice. The extrinsic curvature 
of the Lemaitre coordinate can also be similarly transformed to the induced coordinate. In 
the set of induced coordinates, the momenta P/ r (momentum in the radial direction), P e 7 e 
(momentum in the 9 direction) and P/ (momentum in the direction) were calculated in 
[4]. The two surface bits used to compute the momenta were bits of 2-spheres, and in the 
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limit the area of these bits went to zero, Pi = S e El , where S e is the area of the two 

' e a c a e a ' e a 

pi pi pi 

surface. Given that, qq ab = E al E bl = „ ea * b we get q = det^ = P. Thus 

1 P 1 P 1 

q ab = — — -— ^ — — (14) 

P S ea S £b 

Calculating 

- P 1 P 1 - 

<^|-&-a|^> (15) 

should be enough to calculate corrections to the metric. 

For the specific purpose of calculating corrections to the unstable orbits, we find the 
corrections to the radial metric. 

^MfbH 1 * (i6) 



where P £r = ^jPl r Pl r is the gauge invariant momentum. In the limit the S er — > 0, the 
P Cr = 2r sm ® SdS( l> a nd S er = 2565(f). Needless to say in this approximation, we correctly 
recover q rr = 1 + 0(1). Further, as we consider orbits r > r g , the |/(P er )| [4] gives a 
fractional contribution to the formulas. 

Having found the correction to the induced metric in the r = r c slice, we then per- 
form a coordinate transformation to find the corrections to the metric in the Schwarzschild 
coordinates. The transformations are 

u _ gTT _|_ ^^_ g RR g rr _ ^H^ g rr + dr dr g RR 

dr dr dR dR dr dr dR dR 

r . dr dt dr dt RR . . 

9 rt = -iz-r9 TT + -ns-^9 RR (18) 



We have 



dr dr dR dR^ 



dt 1 dt , /' dr [r~ dr fr g 



dr 1 - /' dR 1 - /' dR V r dr 
This gives, in particular the corrections to g tr as 



-y- = ±<P (19) 



9 tr = (^Y 2 iffe) (20) 



l- Tjl \ r y V^e , 

The reason we are giving this in details is because the quantum gravity effects have 
created a g rt term in the corrected metric which normally wouldn't have been there. The 
g rt of the inverse metric is —g tr /{g tt g rr ) and is thus given by the same rhs of (20). This 
of course diverges at the horizon, but this is a sign of the failure of the coordinates. The 
corrections to the cross terms in the metric g^, g t e, g ri j>, g r 0, 9e^> are not there, as by choice of 
gauge in the internal directions the cross terms like P l e Pi = to order t 
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III. SEMICLASSICAL EFFECTS ON UNSTABLE ORBITS 

In [1], the unstable orbits in the Schwarzschild space-time appear with Lyapunov expo- 
nents. We see how their calculations differ on the inclusion of the quantum gravity correc- 
tions. Initially, we take linearised corrections g M „ + hp V , and then take the as obtained 
in the previous section. We find the most generalised equation of motion, and the type of 
unstable orbits which might emerge from the semiclassical perturbations. 

The geodesic equations are derived in their first order form in phase space. The La- 
grangian for a particle is given by (in the squared form) using the conventions which appear 
in [1] (s is a parameter), 

2 ,, , ( \ / , \ 2 



C 



(21) 



The 9 = 7r/2 in the above. The Equation of motion are derived from the above as: 

: E 



HI 



Ml 



9(r)( d A+h,M 



ds 



ds 



Pr 



(22) 
(23) 



where f(r) = - (l - + h u , g(r) = 
subsequent equations f(r),g(r) = f,g. 



h rr and E is constant of motion. In the 



One can solve (22,23) for dt/ds and dr /ds which has to be used in the equation of motion 

dt E p r h r t 



ds f gf 
dr p r Eh rt 



+ 0(h 2 ) 
+ 0(h 2 ) 



(24) 
(25) 



ds g fg 

The equation of motion for p r is given by = |p and in addition we make a transfer 
mation to the t coordinate by multiplying and keeping terms to 0(h 2 ): 



dp r ds 
ds dt 



drf 



+ d r g(r) 



K ds t 

' dr^ 
ds , 



2d r h tr 

'ds\ 
dt + 



dr' 
ds 



d r q(r) 



ds 



ds 
dt 



(26) 
(27) 
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dp r 

~dt 



'E 
J 



p r h r 

~gj 



+ d r h r t 



P 



2Ep r h H 



IPr 

'E 
J ' 



p r h rt 



p r hrt 



f fg 



28) 



Where = L where q = r 2 sin 2 6 + and L is the angular momentum constant. Using 
the above, one finds the matrix which describes the perturbation flow of the v = {r,p r ) is 

'rom 1,2. Listing 

h rt d r In / d r h rt 



= Kij5vi(t). where i runs from 1,2. Listing the elements of K^, one finds 



K r 



fj 



fj 



Eg 



K 



1 „ „ s h r t d r h r t ld r laqLhrtf 



PrPr 



K 



p r r 



2 g 
l -d 2 .\nfE+ l -d r 



g 2" qgE 2 
d r In qL 2 f \ 



(29) 
(30) 
(31) 
(32) 



In the above we take the r = const orbits or p r = Eh r t/f- The eigenvalues of the matrix are 

1 



^1,2 



E-p r p r ] i \J (K rr ~\~ Kp r pr) ^(E-rrE-p r p r E-rp r K-p r r) 



(33) 



The remarkable fact is that, for orbits which normally have imaginary Lyapunov exponents, 
there is a real part. I thus concentrate on the terms K rr + K PrPr , which might give detectable 
effects. 



K rr + Kp r p r = ^ 



d r hxf , L 2 f 



g 



E 2 qg 



—d r In q 



(34) 



Using the definition for g(r) and f(r) one finds that 



h r t 


'-2- 


L 2 


K r ) r 3 


~2 




-E~ 2 





(35) 

This factor is non-zero, positive and real for all orbits r > r g if we take g r t = h r t for the 
black hole metric. For the white hole metric, the Lyapunov exponent will be positive for 
backwards flow in the time direction. Note, the Lyapunov exponent is postive as the quantum 
fluctuations are negative in the black hole space-time, and the evolution matrix defined as 
Lij(t) = K im L m j(t) has a growing perturbation. The principal Lyapunov exponent is defined 
in terms of the evolution matrix as: A = lim^oo^ In |L^Ljj|. 

Thus for a given fixed value of r, L, E we can find A. If we take a given L, E, the radius of 
the circular orbit in the quantum corrected metric is slightly perturbed. The radius of the 



s 



circular orbit is calculated by finding the minimum of the effective potential which appears 
in the description of the trajectory of an particle. This is obtained in the case of the classical 
metric by using the fact that along a time like trajectory, the distance squared is a number 
(I), 

E 2 1 (drV L 2 

r r \ / 

where in the expression for the line element we have put in the constants E, L by replacing 
the j-, ^ appropriately It follows that this is same as 

(§)' + V<r) = (37) 

where V(r) = (l - r -A + lj + E 2 . The circular orbits are then obtained as a solution 
for f£ = 0. In case of the corrected metric, the same (36) gets modified to the following: 

E 2 (dr\ 2 L 2 



The potential is then identified as 



V h (r) = - 
9 



l 2 e 2 ; 

Q f 



(39) 



Solving for av ^ = at r + S , where r extremises the classical potential 



V(r )h'(r ) - h'(r ) 

s " = (40) 



where ' denotes the derivative wrt r, and h(r) = (1 — h u + -ffi + an d h — (1 + 

L 2 / r 2 )h tt + (l — y-J + E 2 h^/r 2 . However, this correction to the orbit radii does not affect 
the evaluation of the Lyapunov exponent in the linearized approximation. 

Next we plug in the values of r, L, E (33), and find the Lyapunov exponent. We take the 
same circular orbit considered in [1] for the classical Lyapunov exponent in order to see the 
changes in the quantum corrected metric. In units of GM/c 2 as in [1] we use L — 4, E — 1 
and get two orbits at r = 4 and r = 12. In case of r$ = 4, the classical Lyapunov exponent 
gets corrected, but to such a small extent, that it is unlikely to manifest itself. Of interest is 
the r = 12 orbit, which is classically stable. Is there anyway that the quantum instability 
will manifest itself? This is a interesting question and we find that for tq = 12, the Lyapunov 
exponent to be: (The l p = 1.62 x l(T 35 m, G = 6.67 x lO^m^kg^s- 2 , c = 2.99 x 10 8 m/s.) 

A = ik fi b i0 - m7) (41) 
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The t — -|, and one gets after putting in appropriate dimensional units like the speed of 
light 

A = 10 40 " 3n (8.6)/ (^j /a (42) 

The radius r g of the black hole is taken to be 10 n l p . Thus for astrophysical black holes 
where n > 33, this is a number unlikely to manifest itself in any experiments. For black 
holes of horizon radius 10 41 / 3 Z P ~ 10 13 ' 67 / p this is order 1. Thus the chaos time scale for this 
is T\ ~ Is. But in general, the chaos time scale would be 

T a «1.1x10 3 "- 41 T7 ^ tS (43) 

IV. DETECTING THE INSTABILITIES 

We investigate the plausible effects of the positive Lyapunov exponent on physical phe- 
nomena. 



A. Gravity waves and New Instabilities 

The interesting effect of the chaos is in the detection of gravitation waves, as discussed in 
the papers [1] . For primordial black holes, the time scale of chaos is found to be T\ — 1/A 
which is order Is for black holes with horizon radius 10 13 ' 67 / p , and we can make an attempt 
for investigating their effects on gravitational wave forms. The gravitational wave has a 
time scale T w = 27r/0. Given that in the frame of the observer, <f) = L(r — 2)/Er 3 for the 
particular orbit under consideration, we get T u ~ 10 n_42 s which is clearly smaller than the 
chaos time scale for n=13.67. So chaos will not set in. Infact if we take the ratios of the 
time scales we get tP- ~ 10 2n+1 . This is not less than 1 even for n=0 or Planck length black 
holes which are anyway beyond the semiclassical regime. 

We can still try to find the radiative decay rate for a test particle orbiting the black 
hole, and that is given by = 2§6 ^'/T' w h ere = m\rri2{jn\ + m^) where m\ is the mass 
of the black hole and m 2 is the mass of the test particle. For chaos to have any effect on 
the gravitation wave decay rate, we get a bound for the mass of the test particle: m-i < 
8.2 x 10~ n ~ 9 kg. This is obviously a restrictive bound (mass of the electron is ~ I0~ 31 kg) , 
and we are unlikely to find any signatures of the instability in the gravity wave phenomena 
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for astrophysical black holes (for solar mass black holes n > 38). Note that for primordial 
black holes we should be really using quantum corrected formulas for 'gravity wave emission'. 
Further, the Lyapunov exponents for null geodesies can be derived as a limiting case of the 
above calculation, and these determine the quasinormal decay of a black hole [6] . Thus these 
new semiclassical instabilities will create new quasinormal decay modes for primordial black 
holes. 



B. Magnification of quantum gravity corrections to classically unstable orbits. 

As described in the previous section radii of circular orbit get modified by the function 
So (40). We find that the corrections to the unstable orbit radii get magnified in time. To 
see this, we do the following calculation: As per [1], the Lyapunov exponent in appropriate 
units is A ~ 10 43_n / s - Substituting this in the eom we find that, using the principal 

Lyapunov exponent, one has 



Sr(t) = exp (1.09 x 10 42 " n tj S (44) 

Thus a particle traversing in the quantum gravity corrected metric will start deviating 
from its trajectory, and this deviation will grow with time. Let us solve for the time, 

t = 9.17x 10"- 43 ln(^^) S (45) 

Given that S ~ r g t ~ 10 _n / p , for Sr(t) ~ 1 cm, we require: 

t = 9.17 x 10 n ~ 43 ln(10 n+33 )s = 2.1 (n + 33) 10 n ~ 42 s (46) 

Thus, for black holes with horizon radius 10 42 Z p this time is of the order of 1 s. This is in 
the range of astrophysical black holes, and is a very significant effect. Thus any test particle 
traversing in unstable orbits will start deviating from the classical radii spontanously due to 
quantum gravity corrections. However, in an accreting disk the 'instability' due to quantum 
fluctuations will be difficult to isolate from external matter perturbations. Thus though the 
calculations is one of the first attempts to magnify quantum gravity effects to macroscopic 
length scales, the above is still a theoretical result, and an appropriate experiment has to 
be designed to obtain any measurable effect. 
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V. CONCLUSIONS 

In this article we calculated non-zero positive Lyapunov exponents for orbits in a semiclas- 
sically corrected metric of a Schwarzschild space-time. However, the Lyapunov exponents 
are suppressed by the semiclassical parameter which controls the semiclassical corrections. 
The instabilities are almost undetectable for astrophysical black holes, and a particle will 
have stable motion. For black holes which have horizon radius of ~ 10 1A l p the Lyapunov 
exponent is order Is -1 , but there is no detectable effect of these instabilities for gravitational 
waves or radiation decay from particles orbiting the primordial black holes. The instabilities 
might manifest themselves in other astrophysical phenomena and black hole ring down. We 
further make the interesting observation that due to semiclassical corrections to the metric, 
the radii of classical circular orbits get corrected. Normally these corrections would be neg- 
ligible, but the corrections to the classically unstable orbits get magnified in time, and for 
astrophysical black holes the deviations grow to order 1 cm in 1 s! This is a spontaneous 
quantum effect, and not the result of any external perturbations of the black hole metric. 
However a tangible experiment has to be designed to make this spontaneous effect detectable 
in a real accreting situation. Further this effect is based on circular unstable orbits, but the 
parameter space of orbits has to be extended, and study of the behaviour of homoclinic 
orbits which have overlap with the unstable orbits has to be done. 

This calculation can be easily extended to the case of unstable null geodesies, and the 
Lyapunov exponent of null geodesies appears in the calculation of the quasi-normal frequency 
of black hole ring down [6]. Moreover, the Lyapunov exponents for Kerr black holes are 
relevant for black hole mergers [7]. The coherent state has to be defined for the Kerr metric 
in a co-rotating frame, and this is work in progress. We expect to extend the calculation of 
this paper to the case of null geodesies, Kerr black holes and find experiments which will 
detect the effects of quantum gravity corrections to the metric. 
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